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Abstract 

Nanocrystals  of  hexagonal  nickel  carbide  have  been  synthesized  via  physical  vapor  deposition  of 
elemental  nickel  onto  the  surface  of  vertically  aligned  carbon  nanotubes.  Combining  high- 
resolution  transmission  electron  microscopy  (HRTEM)  with  three-dimensional  atom  probe 
tomography  (3DAP)  confirmed  that  these  nanocrystals  have  a  hexagonal  structure,  are  enriched 
in  carbon,  and  have  a  composition  of  ~Ni-25at%C  (M3C).  This  metastable  hexagonal  nickel 
carbide  phase  appears  to  be  stabilized  due  to  the  growth  of  the  nanocrystals  on  the  surface  of  the 
nanotubes  that  act  as  a  template  and  also  as  a  source  of  carbon.  The  stability  of  this  nickel 
carbide  phase  has  also  been  investigated  by  Density  Functional  Theory  (DFT)  calculations  and 
compared  to  the  experimental  results. 

Keywords:  Hybrid  Nanostrucutre;  Carbon  Nanotube;  Metal  Nanocrystal;  Hexagonal  Nickel; 
Atom  probe  tomography. 

Synthesis  of  nanocrystalline  hybrid  materials  based  on  carbon  nanotubes  with  metals  or  oxides 
has  been  of  substantial  interest  in  recent  years  due  to  their  potential  applications  in  bio  sensing, 
energy  storage,  catalysis,  electronic  devices  as  well  as  structural  applications  [1-5].  For  example, 
tin  oxide  (SnCF)  nanocrystals  deposited  by  an  aero-gel  process  onto  multi-walled  carbon 
nanotubes  have  been  found  to  be  promising  for  gas-sensing  applications  at  room  temperatures 
[1]  indicating  that  carbon  nanotubes  can  be  used  as  templates  for  deposition  or  growth  of 
nanocrystals  for  various  applications.  Despite  the  various  reports  on  synthesis  of  carbon 
nanotubes  coated  with  nanocrystals,  there  have  been  rather  limited  high-resolution  investigations 
of  the  structural  and  compositional  details  associated  with  these  nanocrystals.  Similarly  there  has 
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been  a  lot  of  interest  in  the  synthesis  of  nanocrystalline  metallic  materials  such  as  nickel,  due  to 
their  application  in  a  wide  range  of  engineering  fields,  specifically  catalysis,  energy  absorption, 
and  magnetic  recording  media  [6-8].  Elemental  nickel,  in  its  equilibrium  form  at  room 
temperature  exhibits  a  face  centered  cubic  (FCC)  structure.  Interestingly,  there  have  been  reports 
in  the  literature  related  to  the  stabilization  of  a  non-equilibrium  hexagonal  structure  (hexagonal 
close-packed  or  HCP)  in  nanostructured  forms  of  nickel,  including  HCP  metastable  nickel 
carbide  nanocrystals  [9-12].  The  HCP  phase  is  typically  observed  in  case  of  far-from  equilibrium 
synthesis  routes  such  as  wet  chemical  synthesis  at  very  low  temperatures  when  the  nanocrystals 
grow  from  an  amorphous  phase  [9],  decomposition  of  neutral  metal  complex  [10],  thin  film 
processing  via  hetero  epitaxial  growth  [11],  and  mechanical  alloying  via  high  energy  ball  milling 
[12].  These  synthesis  methods  indicate  that  HCP  nickel  is  only  observed  below  a  certain  critical 
crystallite  size  (~  4  nm)  and  changes  to  the  equilibrium  FCC  structure  for  larger  crystallite  sizes 
[9-10].  Nickel  does  not  form  a  carbide  under  equilibrium  conditions,  however,  there  have  been 
reports  of  metastable  HCP  nickel  carbides  (NFC)  forming  as  a  result  of  mechanical  alloying  in 
high  energy  ball-milled  mixtures  [12].  These  HCP  nickel  carbide  (NFC)  nanocrystals  also 
exhibited  a  highly  disordered  and  defective  internal  structure  due  to  extreme  local  pressures  and 
temperatures  involved  in  far-from  equilibrium  processing,  such  as  high  energy  ball  milling  [12]. 
Furthermore,  nanocrystals  synthesized  from  a  carbon-containing  neutral  metal  complex  of 
nickel,  indicate  that  the  HCP  phase  can  be  stabilized  in  much  larger  nanocrystals  (~  20  nm), 
possibly  due  to  the  presence  of  carbon  [10].  However,  a  detailed  investigation  and  clear 
understanding  of  the  role  of  carbon  on  the  stabilization  of  the  HCP  phase  in  nickel  nanocrystals 
is  presently  lacking. 

This  paper  reports  the  templated  synthesis  of  hexagonal  nickel  carbide  nanocrystals  via  physical 
vapor  deposition  (PVD)  of  elemental  nickel  onto  the  surface  of  double- walled  carbon  nanotubes. 
PVD  processes  such  as  sputter  deposition  offer  a  relatively  simple  and  clean  route  for 
synthesizing  nanocrystals  on  nanotube  surfaces.  Electron  diffraction  and  high  resolution 
transmission  electron  microscopy  (HRTEM)  was  used  to  identify  the  crystal  structure  and 
morphologies,  while  three-dimensional  atom  probe  tomography  was  used  to  determine  the 
composition  of  the  nickel  carbide  nanocrystals.  The  coupling  of  these  two  advanced 
characterization  techniques  permitted  the  determination  of  both  the  atomic  scale  structure  and 
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composition  of  these  nanocrystals  for  the  first  time.  Furthermore,  synergistic  coupling  of  these 
experimental  results  with  first-principles  electronic  structure  based  calculations  of  the  relative 
energies  of  different  metastable  nickel  and  nickel  carbide  phases  has  led  to  a  better  insight  into 
the  role  of  carbon  and  the  template  offered  by  the  nanotube  surface  on  the  formation  of 
metastable  nickel  carbide  nanocrystals. 

2.  Experiment 

Using  sputtering,  elemental  nickel  was  deposited  onto  vertically  aligned  CNT  in  a  PVD  system. 
Pure  elemental  Ni  (99.99%)  targets  from  Kurt  Lesker  Co.,  were  used  for  the  deposition.  Si 
wafers  with  a  high  density  of  vertically  grown  CNTs  on  the  surface  were  used  as  the  substrate 
for  the  Ni  depostion.  A  10  mm  x  10  mm  piece  of  Si  wafer,  coated  with  the  CNTs,  was  loaded 
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into  the  PVD  chamber  and  pumped  down  to  a  base  pressure  ~  10"  torr.  Pure  Argon  sputtering 
gas  was  introduced  into  the  chamber  to  achieve  a  working  pressure  ~  5  mTorr.  Ni  was  deposited 
using  the  DC  magnetron  sputtering  power  at  200W  at  room  temperature  for  60  minutes.  In  order 
to  characterize  the  nanostructured  Ni  deposited  on  the  CNT  surface,  site-specific  samples  were 
prepared  for  both  transmission  electron  microscopy  (TEM)  and  3D  atom  probe  (3DAP)  studies, 
using  the  FEI  Nova  200  NanoLab  -  dual  beam  FIB  instrument.  In  order  to  minimize  the  Gallium 
ion  damage  during  the  sample  preparation,  the  area  of  interest  was  protected  by  depositing  a  Pt 
overlayer  on  the  region  of  interest  and  subsequently  lifting  out  samples.  Lifted-out  samples  were 
attached  to  a  copper  grid  for  TEM  studies.  Additional  thinning  and  cleaning  to  obatin  less  than 
50  nm  thickness  using  the  FIB  was  done  at  30  KeV  and  5  KeV  consequently  to  remove  the 
redeposition  and  ion  beam  damage.  The  site  specific  samples  have  been  characterized  using  the 
FEI  Tecnai  F20  field  emission  gun  TEM  operating  at  200  kV.  The  sample  for  atom  probe 
analysis  was  also  prepared  using  the  dual-beam  FIB.  3DAP  experiments  were  conducted  on  a 
LEAP  3000  Local  Electrode  Atom  Probe  (LEAP™)  microscope  system  manufactured  by  Imago 
Scientific  Instruments,  Inc.  The  laser-pulse  evaporation  method  was  used  for  analyzing  these 
samples  at  a  temperature  of  60K  with  0.5%  average  evaporation  rate  and  a  laser  pulse  energy  of 
0.5  nJ. 
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3.  Results  and  Discussion 


Fig.  1(a)  shows  a  scanning  electron  microscopy  (SEM)  image  on  the  as-received  vertically 
aligned  carbon  nanotubes  grown  by  chemical  vapor  deposition  method  on  a  Si  wafer.  The  high 
density  of  aligned  CNTs  have  a  double-walled  structure  with  a  diameter  ~  4  nm  and  length  ~  650 
pm  length.  The  detailed  properties  of  these  nanotubes  are  reported  elsewhere  [13].  Fig.  1(b) 
shows  a  SEM  image  from  the  nickel-coated  carbon  nanotubes,  clearly  exhibiting  the  granular 
morphology  of  nanocrystals  uniformly  coating  the  surface  of  the  nanotubes  (higher 
magnification  shown  as  an  inset).  From  bright-field  TEM  images,  such  as  the  representative  one 
shown  in  Fig.  1(c),  the  average  crystalline  size  was  determined  to  be  ~5  nm.  The  electron 
diffraction,  comprising  multiple  rings,  shown  as  an  inset  in  Fig.  1(c),  confirmed  the 
polycrystalline  nature  of  the  vapor  deposited  nanocrystals  on  the  surface  of  the  nanotubes. 

For  a  detail  analysis  of  morphology  and  crystal  structure  of  these  nanocrystals,  an  HRTEM 
image  as  well  as  an  enlarged  selected  area  diffraction  (SAD)  pattern  are  shown  in  Figs.  2(a)  and 
(b)  respectively.  Fig.  2(a)  clearly  shows  the  formation  of  nanocrystals  of  relatively  uniform  size 
on  the  surface  of  the  nanotubes,  with  the  underlying  graphene  layers  visible  in  some  cases.  The 
nanocrystals  were  typically  less  than  ~  3nm  in  size.  Most  of  the  primary  rings  observed  in  the 
SAD  pattern  shown  in  Fig.  2(b)  can  be  consistently  indexed  based  on  a  hexagonal  (HCP) 
structure  with  a  single  weak  ring  corresponding  to  the  face-centered  cubic  (FCC)  structure  of  Ni. 

o 

The  lattice  parameter  of  the  HCP  phase  was  determined  to  be  a  =  2.6  ,  c  =4.3A  and  found  to  be 
in  close  agreement  with  those  reported  previously  by  Tian  et.  al  for  hexagonal  nickel 
nanocrystals  [11]  as  well  as  by  Yue  et.  al.  for  metastable  M3C  nanocrystals  [12].  The  equilibrium 
phase  diagram  does  not  exhibit  either  any  hexagonal  Ni  or  nickel  carbide  phase.  It  appears  that 
the  carbon  nanotubes  can  affect  the  formation  of  this  hexagonal  nickel  phase  by  providing  a 
template  for  heterogeneous  nucleation  during  the  physical  vapor  depostion  process.  Fig.  2(c) 
shows  a  relatively  large  Ni  nanocrystal  (-20  nm)  on  the  CNT  surface  exhibiting  a  hexagonal 

structure  as  revealed  by  the  corresponding  [  1213  ]  zone  axis  SAD  pattern  (inset  in  Fig.  2(c)).  The 
evidence  of  a  hexagonal  structure  in  this  relatively  large  nanocrystal  contradicts  the  critical  size 
(~  4nm)  of  FCC-HCP  transition  reported  previously  for  Ni  nanocrystals  by  Illy  et  al.  [9]. 
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Figs.  3(a)  and  3(b)  are  SEM  images,  respectively  showing  the  lift-out  procedure  and  the  sample 
prepared  for  3DAP  studies  of  the  Ni  nanocrystals  deposited  on  the  carbon  nanotubes.  The  lifted- 
out  sample  was  attached  on  a  Si  micropillar,  and  then  subsequently  milled  using  the  Ga  ion  beam 
with  a  final  tip  radius  less  than  30nm  as  shown  in  Fig.  3(b).  The  3DAP  sample  preparation  using 
FIB  was  well  described  elsewhere  [14].  Fig.  3(c)  shows  a  3D  atom  probe  reconstruction  of  the 
sample  shown  in  Fig.  3(b).  The  C  atoms  have  been  plotted  as  red  (color  online)  dots  along  with  a 
solid  isoconcentration  surface  (isosurface  in  short)  in  blue  (color  online)  for  Ni=50at%.  Even 
though  the  Ni  isosurface  shows  a  clear  interface  between  the  Ni  nanocrystal  and  CNT  sides,  the 
composition  profiles  across  the  nickel  and  CNT  interface,  plotted  in  Fig.  3(d),  exhibit  substantial 
intermixing  between  the  two  regions.  The  average  composition  of  carbon  in  the  Ni  coated  area 
(nickel  nanocystals)  is  -25  at.%,  which  corresponds  to  ~Ni3C  stoichiometry  (refer  Fig.  3(d)). 
Interestingly,  while  the  metastable  Ni3C  phase  has  been  proposed  to  form  in  high-energy  ball- 
milled  nanocrystalline  nickel  samples,  there  was  no  compositional  evidence  presented  in  support 
of  this  argument  [12].  However,  the  reported  HCP  crystal  structure  and  lattice  parameters  for  this 
metastable  carbide  are  in  good  agreement  with  the  results  of  the  present  study.  Yue  et.  al 
suggested  carbon  seems  to  be  a  significant  role  on  the  stabilization  of  HCP  Ni  phase,  especially 
at  lower  growth  temperatures.  The  present  study  suggest  that  hexagonal  Ni3C  nanocrystals  can 
be  stabilized  up  to  sizes  as  large  as  ~20nm  when  grown  on  the  surface  of  carbon  nanotubes. 
Therefore,  it  appears  that  the  nanotubes  not  only  provide  a  template  for  the  growth  of  the 
hexagonal  nickel  phase  but  also  increase  the  stability  of  this  phase  due  to  the  diffusion  of  carbon 
from  the  nanotubes  into  the  nanocrystals.  Under  these  circumstances  it  becomes  imperative  to 
investigate  the  relative  stabilities  of  HCP  Ni  and  HCP  Ni3C  phases  with  respect  to  the 
equilibrium  FCC  Ni  phase.  Such  an  analysis  has  been  carried  out  using  density  functional  theory 
(DFT)  as  described  below. 

We  have  performed  plane  wave  Density  Functional  Theory  (DFT)  total  energy  calculations  of 
HCP  and  FCC  Ni,  as  well  as  orthorombic  (cementite)  and  HCP  Ni3C  [16].  Equation  of  states  of 
these  structures  were  obtained  by  calculating  the  total  energy  versus  volume  curves  with  unit  cell 
volumes  spanning  10%  up  and  below  those  of  the  experimental  structures.  The  cell  shape  and 
atom  positions  were  fully  relaxed  within  the  theory  of  DFT  at  each  of  the  volumes.  The  final 
equation  of  states  curves  are  shown  in  Fig.  4.  Generalized  gradient  exchange  correlation 
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functionals  with  the  PBE  form  [17]  and  Projected  Augmented  Wave  (PAW)  [18] 
pseudopotentials  to  describe  the  electron  -  ion  core  interactions  were  used  in  the  DFT 
calculations  as  implemented  in  the  Vienna  ab  initio  Simulation  Package  (VASP)  [19-20].  Ni  3s 
and  3p  electrons  were  treated  as  valence  electrons  (3s  p  d  4s“).  Plane  wave  kinetic  energy  cutoff 
of  400  eV  and  9x9x9  for  NOC  and  21x21x21  for  Ni  K-point  sampling  of  Brillouin  zone  were 
used  in  the  calculations.  Geometry  optimizations  were  performed  using  the  congugaed  gradient 
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method  until  forces  acting  on  each  atom  being  smaller  than  0.01  eV/A.  Spin  polarized 
calculations  was  employed  in  all  geometry  optimzation  and  energy  calculations  due  to  the 
maganetic  nature  of  Ni.  HCP  M3C  was  found  to  have  zero  magnetization  but  both  HCP  and  FCC 
Ni  exhibited  ferromagnetic  behavior  with  residual  magnetic  moments  of  0.63  pB- 

DFT  calculations  show  that  the  FCC  Ni  is  more  stable  than  HCP  Ni  by  0.024  eV/atom  at 
equilibrium  geometry,  consistent  with  experiments  and  previous  DFT  calculations  and  verifies 
that  spin  polarized  calculations  are  important  in  nickel  and  related  systems  [22].  In  Fig.  4,  in 
order  to  be  comparable  to  the  M3C  equation  of  state,  the  energies  and  volumes  of  HCP  and  FCC 
Ni  are  ploted  for  three  Ni  atoms  and  the  energy  was  shifted  by  adding  the  energy  per  atom  of 
carbon  (in  the  diamond  structure)  from  calculations  of  the  same  conditions  (-9.110  eV/atom). 
The  orthorhombic  cementite  structured  M3C  equation  of  state  was  also  calculated  (not  shown) 
and  it  is  found  that  it  has  a  higher  energy  than  the  HCP  structure  by  0.157  eV  per  M3C  unit 
which  is  in  good  agreement  with  experimental  observations  of  the  existence  of  only  hexagonal 
Ni3C.  The  energy  versus  volume  curves  were  fitted  to  Birch-Mumaghan  equation  of  state  and 
equilibrium  volume,  cell  parameter  and  bulk  modulus  were  obtained.  For  hexagonal  M3C,  the 
equilibrium  volume  obtained  from  equation  of  state  is  about  2.82%  larger  than  the  experimental 
volume  determined  using  electron  diffraction  and  the  corresponding  unit  cell  sizes  are  0.93% 
larger.  The  bulk  modulus  from  equation  of  state  is  188.2  GPa  for  NBC.  The  equilibrium  M3C 
structure  was  used  to  calculate  the  electron  diffraction  patterns.  The  calculated  TEM  pattern 
agrees  well  with  experimental  data  both  in  terms  of  the  general  hexagonal  symmetry  and  the 
interatomic  spacing  of  the  planes  in  the  diffraction  pattern  as  shown  in  Table  1. 

The  formation  of  hexagonal  M3C  can  become  favorable  in  the  consideration  of  total  energy  of 
the  Ni/CNT  systems  based  on  our  DFT  calculations.  The  energy  difference  between  the  HCP  Ni 
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+  carbon  and  the  hexagonal  N13C  is  0.063  cV/Ni^C  unit.  Here  the  energy  of  carbon  used  for 
adding  to  the  HCP  Ni  value  is  for  the  diamond  structure,  which  compared  to  graphite  is  lower  by 
0.0197  eV/atom.  However,  due  to  the  fact  that  in  this  case  the  carbon  atoms  are  in  the  form  of  a 
nanotube,  the  total  energy  of  carbon  can  be  considerably  higher  than  that  of  diamond  or  graphite. 
It  was  found  from  DFT  calculations  that  the  energy  of  formation  of  carbon  nanotubes  increases 
with  decreasing  tube  diameters  [23].  This  will  shift  the  equation  of  state  for  FCC  and  HCP  Ni  + 
C  to  a  higher  energy  level,  which  will  eventually  favor  the  formation  of  HCP  N  i  3C  nanocrystals 
rather  than  HCP  Ni  nanocrystals  on  the  surface  of  the  CNTs.  Furthermore,  interfacial  energy 
considerations  are  also  likely  to  aid  in  the  stabilization  of  the  HCP  NFC  nanocrystals. 

Sputter  deposition  of  elemental  Ni  onto  vertically-aligned  carbon  nanotubes  results  in  the 
templated  nucleation  and  growth  of  NFC  nanocrystals  exhibiting  a  hexagonal  crystal  structure. 
The  nanotubes  not  only  act  as  a  template  for  the  synthesis  of  these  nanocrystals  but  also  act  as  a 
source  of  carbon.  The  hexagonal  phase  can  be  stabilized  in  these  nanocrystals  to  relatively  large 
sizes  (~  20  nm).  Density  functional  theory  calculations  indicate  that  the  HCP  N i 3C  phase  is 
metastable  with  respect  to  FCC  Ni  +  C  (in  diamond).  However,  the  total  energy  difference 
between  HCP  NFC  and  HCP  Ni  +  C  (in  diamond)  is  relatively  small  and  it  is  quite  likely  that  the 
HCP  NFC  phase  might  have  a  lower  energy  than  HCP  Ni  +  C  (in  nanotubes),  resulting  in  their 
stabilization  on  the  nano  tube  surface. 

Acknowledgements 

This  work  was  financially  supported  by  U.  S.  Air  Force  Research  Laboratory  (AFRL  ISES 
contract)  and  the  U.  S.  Air  Force  Office  of  Scientific  Research  (AFOSR  Grant  #  FA9550-06-1- 
0193).  The  authors  also  gratefully  acknowledge  the  Center  for  Advanced  Research  and 
Technology  (CART)  at  the  University  of  North  Texas. 


7 


[1]  G.  Lu,  L.E.  Ocola,  and  J.  Chen,  Adv.  Mater.  21,  1  (2009). 

[2]  D.S.  Kim,  S.M.  Lee,  R.  Scholz,  M.  Knez,  U.  Gosele,  J.  Fallert,  H.  Kalt,  and  M.  Zacharias, 
Appl.  Phys.  Lett.  93,  103108  (2008). 

[3]  L.  Zhu,  G.  Lu,  and,  J.  Chen,  J.  Heat  Transfer,  130,  044502-1  (2008). 

[4]  EJ.  Garcia,  A.J.  Hart,  B.L.  Wardle,  and  A.H.  Slocum,  Adv.  Mater.  19,  2151  (2007). 

[5]  K.  Ota,  A.  Kawabata,  H.  Murakami,  E.  Kita,  Mater.  Trans,  42,  1684  (2001). 

[6]  N.  Cordente,  M.  Respaud,  F.  Senocq,  M.J.  Casanove,  C.  Amiens,  and  B.  Chaudret,  Nano 
Lett.  1,  565  (2001). 

[7]  J.C.  Bao,  C.Y.  Tie,  Z.  Xu,  Q.F.  Zhou,  D.  Shen,  Q.  Ma,  Adv.  Mater.  13.  1631  (2001). 

[8]  Z.P.  Liu,  S.  Li,  Y.  Yang,  S.  Peng,  Z.K.  Hu,  Y.T.  Qian,  Adv.  Mater.  15.  1946  (2003). 

[9]  S.  Illy,  O.  Tillement,  F.  Machizud,  J.M.  Dubois,  F.  Massicot,  Y.  Fort,  and  J.  Ghanbaja,  Phil. 
Mag.  A,  79,  1021  (1999). 

[10]  V.  Rodriguez-Gonzalez,  E.  Marceau,  P.  Beaunier,  M.  Che,  C.  Train,  J.  Solid  State  Chem. 
180,  22  (2007). 

[11]  W.  Tian,  H.P.  Sun,  X.Q.  Pan,  J.H.  Yu,  M.  Yeadon,  C.B.  Boothroyd,  Y.P.  Feng,  R.A. 
Lukaszew,  and  R.  Clarke.  Appl.  Phys.  Lett.  86,  131915  (2005). 

[12]  L.  Yue,  R.  Sabiryanov,  E.M.  Kirkpatrick,  D.L.  Leslie-Pelecky,  Phys.  Rev.  B.  62,  8969 

(2000). 

[13]  X.  Zhang,  Q.  Li,  Y.  Tu,  Y.  Li,  J.  Y.  Coulter,  L.  Zheng,  Y.  Zhao,  Q.  Jia,  D.E.  Peterson,  Y. 
Zhu,  Small,  3.  244  (2007). 

[14]  K.  Thompson,  D.  Lawrence,  D.J.  Larson,  J.D.  Olson,  T.F.  Kelly,  and  B.  Gorman, 
Ultramicroscopy,  107,  131  (2007). 

[15]  M.K.  Miller,  K.  R.  Russell,  K.  Thompson,  R.  Alvis,  and  D.J.  Larson,  Microsc.  Microanal. 
13,  428  (2007). 

[16]  S.  Nagakura,  J.  Phys.  Soc.  Jpn.  13,  1005  (1958). 

[17]  G.  Kresse  and  J.  Hafner,  Phys.  Rev.  B,  49,  14251  (1994). 

[18]  G.  Kresse  and  J.  Furthmuller,  Phys.  Rev.  B,  54,  11169  (1996). 

[19]  J.  P.  Perdew,  K.  Burke,  and  M.  Emzerhof,  Phys.  Rev.  Lett.,  77,  3865  (1996). 

[20]  P.  E.  Blochl,  Phys.  Rev.  B,  50,  17953  (1994). 

[21]  G.  Kresse  and  D.  Joubert,  Phys.  Rev.  B,  59,  1758  (1999). 

[22]  G.  Y.  Guo  and  H.  H.  Wang,  Chin.  J.  Physics,  38,  949  (2000). 

[23]  M.F.  Budyka,  T.S.  Zyubina,  A.G.  Ryabenko,  S.H.  Lin,  A.M.  Mebel,  Chem.  Phys.  Lett.  407, 
266,(2005). 


8 


50  urn 


20  nm 


Fig.  2. 


10 


Fig.  3. 


11 


-25.6 


-25.8 


-26.2 


-26.4 


-26.6 


30  35  40  45 

Volume  (A3) 


Fig.  4  Energy  versus  volume  curves  of  HCP  Ni,  FCC  Ni  and  HCP  NprC  calculated  using  spin 
polarized  periodic  planewave  DFT  with  GGA-PBE  functional.  Open  circle:  N13C;  filled  circle: 
HCP  Ni,  filled  square  FCC  Ni.  For  HCP  and  FCC  Ni,  the  curves  were  plotted  for  3  Ni  atoms  and 
the  energy  shifted  by  adding  the  energy  of  carbon  in  the  diamond  form  to  be  comparable  with 
Ni3C.  Inset:  the  unit  cell  of  FCC  and  HCP  Ni,  and  HCP  N^C.  Grey  ball:  carbon,  blue  ball: 
nickel. 


Table  1.  Comparison  of  the  unit  cell  parameter  of  hexagonal  Ni3C  and  hexagonal  Ni. 


Cell  parameter 

Angles 

Space 

group 

Interlayer 

distances 

TEM 

measured 

distance 

Nl3C  (HCP) 

4.600/13.016 

90.0/90.0/120.0 

R-3C 

10-14:  2.52 

2.5 

( 4.553/12.920 f 

( 90.0/90.0/120.0 )a 

(167) 

11-23:  2.03 

2.0 

11-26:  1.57 

1.5 

Ni  (HCP) 

2.488/4.099 

90.0/90.0/120.0 

P63/MMC 

(2.48/4.08 f 

(90.0/90.0/120.0) 

(194) 

a  experimental  structure  from  ref.  16. 
b  equilibrium  strucure  from  DFT  calculations  from  ref.  22. 
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